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Abstract 

We present in analytic form the matching conditions for the strong coupling 
constant aj lf (/i) at the flavour thresholds to three loops in the modified minimal- 
subtraction scheme. Taking into account the recently calculated coefficient fa of 
the Callan-Symanzik beta function of quantum chrornodynamics, we thus derive a 
four-loop formula for cTC/O together with appropriate relationships between the 
asymptotic scale parameters A( ra -f) for different numbers of flavours rif. 
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The strong coupling constant ai nf \fi) = gl/{ 47t), where g s is the gauge coupling of 
quantum chromodynamics (QCD), is a fundamental parameter of the standard model of 
elementary particle physics; its value af\Mz) is listed among the constants of nature in 
the Review of Particle Physics [Q. Here, /i is the renormalization scale, and rif is the 

number of active quark flavours q, with mass rn q -C fi. The fi dependence of oiP ! \n) is 
controlled by the Callan-Symanzik beta function of QCD, 
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The calculation of the one-loop coefficient f3j } nf> about 25 years ago 0 has led to the 
discovery of asymptotic freedom and to the establishment of QCD as the theory of strong 
interactions. In the class of schemes where the beta function is mass independent, which 
includes the minimal subtraction (MS) schemes of dimensional regularization 0, f3 ’ Q n/ 
and (3^1^ [0 are universal. The results for (3^^ 0 and (3^^ [|] are available in the 
modified MS (MS) scheme 0. For the reader’s convenience, (3^^ (N = 0, ...,3) are 
listed for the rif values of practical interest in Table [T]. 
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In MS-like renormalization schemes, the Appelquist-Carazzone decoupling theorem || 
does not in general apply to quantities that do not represent physical observables, such as 
beta functions or coupling constants, he., quarks with mass m q /i do not automatically 
decouple. The standard procedure to circumvent this problem is to render decoupling 
explicit by using the language of effective field theory. As an idealized situation, consider 
QCD with rii = rif— 1 massless quark flavours and one heavy flavour h, with mass rrih H- 
Then, one constructs an effective rq-fiavour theory by requiring consistency with the full 
n j-fiavour theory at the heavy-quark threshold = 0(rrih). This leads to a nontrivial 
matching condition between the couplings of the two theories. Although, oiy l \mh) = 
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ois nf \rrih) at leading and next-to-leading order, this relation does not generally hold at 
higher orders in the MS scheme. If the // evolution of aP f \/i) is to be performed at N + 1 
loops, he., with the highest coefficient in Eq. ([!]) being P^ f \ then consistency requires 
that the matching conditions be implemented in terms of IV-loop formulae. Then, the 
residual /i dependence of physical observables will be of order IV+ 2 . A pedagogical review 
of the QCD matching conditions at thresholds to two loops may be found in Ref. [|9]. 

The literature contains two conflicting results on the two-loop matching condition for 
a^ lf \fj) in the MS scheme [|TU], [TT[] . The purpose of this letter is to settle this issue by 
an independent calculation and to take the next step, to three loops. As a consequence, 
Eq. (9.7) in the encyclopedia by the Particle Data Group |TJ will be corrected and extended 
by one order. We shall also add the four-loop term in the formula (9.5a) for in 

Ref. 0. 

In order to simplify the notation, we introduce the couplant <A n/ )(/i) = aP t \n)/ 7 T and 
omit the labels // and rtf wherever confusion is impossible. Integrating Eq. (0) leads to 
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where b n = Pn/Po (AI = 1,2,3), A is the so-called asymptotic scale parameter, and C 
is an arbitrary constant. The second equation of Eq. (j2|) is obtained by expanding the 
integrand. The conventional MS definition of A, which we shall adopt in the following, 


Iteratively solving Eq. ( 0 ) yields 
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where L = ln(/i 2 /A 2 ) and terms of 0(1/L 5 ) have been neglected. Equation @ extends 
Eq. (9.5a) of Ref. |T[ to four loops. 

The particular choice of C [0, [12|] in Eq. (|) is predicated on the grounds that it 
suppresses the appearance of a term proportional to (const. /L 2 ) in Eq. (0). For practical 
applications, it might be more useful to define C by equating the one- and two-loop 
expressions of a s (/a), he., by nullifying the 0(1/L 2 ) term in Eq. (0), at some convenient 
reference scale )jl 0 JR|, e.g., at /iq = Mz■ By contrast, in the standard approach, one 


has /i(j = y^eA, which is in the nonperturbative regime. This would lead to the choice 
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C = [bi/(5 q) ln[/3 0 1 ii(/Uq/A 2 )]. The advantage of this convention would be that the values 
of A would be considerably more stable under the inclusion of higher-order corrections. 
Another interesting alternative is to adjust C in such a way that A becomes nf independent 
U . 

ft is interesting to quantitatively investigate the impact of the higher-order terms of 
the beta function in Eq. © on the /i dependence of for fixed nj. For illustration, we 
consider, as an extreme case, the evolution of a[ 5 \/d) from fi = M z down to scales of the 
order of the proton mass. Specifically, we employ the four-loop formula (||) and its A r -loop 
approximations, with AT = 1, 2, 3, which emerge from Eq. ([|) by discarding the terms of 
0(1/L N+1 ). In each case, we determine A® from the condition that af\M z ) = 0.118 [IJ 
be exactly satisfied. For comparison, we also consider the exact solution of Eq. (jl j with all 
known beta-function coefficients included. In Fig. |l]. the various results for 1 / af\ii) are 
plotted versus [i/M z using a logarithmic scale on the abscissa. Consequently, the one-loop 
result appears as a straight line. All curves precisely cross at /i = M z , outside the figure. 
We observe that, for N increasing, the expanded Wloop results of Eq. ([!]) gradually 
approach the exact four-loop solution of Eq. ([[]) in an alternating manner. Down to 
rather low scales, the two-loop result already provides a remarkably useful approximation 
to the exact four-loop result, while the one-loop result is far off. 

Next, we outline the derivation of the three-loop matching condition. In the following, 
unprimed quantities refer to the full n/-flavour theory, while primed objects belong to the 
effective theory with m = rif — 1 flavours. Futhermore, bare quantities are labelled by 
the superscript 0. We wish to derive the decoupling constant ( g in the relation g' s = ( g g s 
between the renormalized couplings g s and g' s . Exploiting knowledge 0 of the coupling 
renormalization constant Z g within either theory, this task is reduced to finding = 

QZ’JZ,. The Ward identity = C?/ (Cs +F) • wllere 


G? = \OK, c® 



A" = CfSv/CSA", 
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with G afi , c, and A^ being the fields of the gluon and the Faddeev-Popov ghost, and the 
Gcc vertex, respectively, then leads us to consider the heavy-quark contributions to the 
corresponding vacuum polarizations and vertex correction, 11^(^), n^(g^), and T g (q c , qc)- 
Specifically, we have 
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In total, we need to compute 3 + 1 + 5 two-loop and 189 + 25 + 228 three-loop Feynman 
diagrams. The 5 two-loop diagrams pertinent to Ci add up to zero. Typical three-loop 
specimen are depicted in Fig. [2|. In order to cope with the enormous complexity of the 
problem at hand, we make successive use of powerful symbolic manipulation programs. 
We generate and compute the relevant diagrams with the packages QGRAF JT5J and 
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Figure 1: /i dependence of a^ s 5 \fi) calculated from af\Mz) = 0.118 using Eq. (|j) at one 
(coarsely dotted), two (dashed), three (dot-dashed), and four (solid) loops. The densely 
dotted line represents the exact solution of Eq. m at four loops. 


MATAD [|16j, respectively. The cancellation of the ultraviolet singularities, the gauge- 
parameter independence, and the renormalization-group (RG) invariance serve as strong 
checks for our calculation. 

If we measure the matching scale in units of the RG-invariant MS mass fih = 
m^/Pi), our result for the ratio of a' = to a = ^"^(/i^) reads 
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Figure 2: Typical three-loop diagrams pertinent to n(t(0;)j U^(q 2 ), and T^(g c , q 5 ). Loopy, 
dashed, and solid lines represent gluons G, Faddeev-Popov ghosts c, and heavy quarks /r, 
respectively. 


where £ h = ln[(0 n/ ^) 2 //4] and 
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Here, ( is Riemann’s zeta function, with values £(2) = 7 t 2 /6 and ((3) ~ 1.202 057. Our 
result for c 2 agrees with Ref. |]TT[, while it disagrees with Ref. 0. For the convenience of 
those readers who prefer to deal with the pole mass Mh, we list here a simple formula [|17| 
for Hh in terms of Mh and A = a^ nf \Mh), which incorporates the well-known two-loop 
relation between mh(M h ) and M h |L8[|. It reads 
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Using a similar relation, with A expressed in terms of a and £/, = In[(// n ^) 2 /M 2 ], we 
may rewrite Eq. (^|) as 
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Going to higher orders, one expects, on general grounds, that the relation between 
o 4 n *V) an d where //' -C /d nf ' > -C /x, becomes insensitive to the choice of /d n /) 

as long as //= 0(rrih)- This has been checked in Ref. 0 for three-loop evolution in 
connection with two-loop matching. Armed with our new results, we are in a position to 
explore the situation at the next order. As an example, we consider the crossing of the 
bottom-quark threshold. In particular, we wish to study how the /b 5 ) dependence of the 
relation between a^ 4 \M T ) and af\M z ) is reduced as we implement four-loop evolution 
with three-loop matching. Our procedure is as follows. We first calculate (/j,^ 5 ) with 
Eq. (|3|) by imposing the condition a^(M T ) = 0.36 ||, then obtain from Eq. (|9|), 

and finally compute af\M z ) with Eq. For consistency, N-loop evolution must be 
accompanied by (N — l)-loop matching, he., if we omit terms of 0(1/L N+1 ) in Eq. (j3[), 
we need to discard those of 0(a N ) in Eq. (||) at the same time. In Fig. [|. the variation of 
af\M z ) with fi^/M b is displayed for the various levels of accuracy, ranging from one- 
loop to four-loop evolution. For illustration, /x® is varied rather extremely, by almost 
two orders of magnitude. While the leading-order result exhibits a strong logarithmic 
behaviour, the analysis is gradually getting more stable as we go to higher orders. The 
four-loop curve is almost flat. Besides the ^ dependence of af\M z ), also its absolute 
normalization is significantly affected by the higher orders. At the central scale /d 5 ) = M b , 
we again encounter an alternating convergence behaviour. We notice that the four-loop 
result is appreciably smaller than the three-loop result, by almost 0.001. This difference 
is comparable in size to the shift in the value of af\M z ) extracted from the measured 
Z-boson hadronic decay width due to the inclusion of the known three-loop correction to 
this observable |T9| . 

As we have learned from Fig. || : in higher orders, the actual value of /d n d does not 
matter as long as it is comparable to the heavy-quark mass. In the context of Eq. (|J), the 
choice /d n d = l^h 0 is particularly convenient, since it eliminates the RG logarithm i b - 
With this convention, we obtain from Eqs. (2j), (|j), and ([]) a simple relationship between 
A' = AGO and A = A("/), viz 
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where lh = ln(p|/A 2 ). The 0(1/1%) term of Eq. (PH) represents a new result. Leaving 
aside this term, Eq. ([□]) disagrees with Eq. (9.7) of Ref. [|TJ. This disagreement may 
partly be traced to the fact that the latter equation is written with the c 2 value obtained 
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a s (M z ) 



Figure 3: /i® dependence of af\M z ) calculated from a^\M T ) = 0.36 and M b = 4.7 GeV 
using Eq. (|3|) at one (dotted), two (dashed), three (dot-dashed), and four (solid) loops in 
connection with Eq. (9j) at the respective order. 


in Ref. Jl0[], which differs from the value listed in Eq. (|7[). Furthermore, in the same 
equation, the terms involving /? 2 should be divided by 4. Equation ([Tlj) represents a closed 
three-loop formula for in terms of and For consistency, it should be used 
in connection with the four-loop expression (3j) for ai n/ ^( / u) with the understanding that 
the underlying flavour thresholds are fixed at = Hh- The inverse relation that gives 
AG/) as a function of A^ n ^ and /i/,. emerges from Eq. ([□]) via the substitutions A <->■ A'; 
Pn (3'n f° r N — 0; • • • j 3; and c n —> ~Cn for N = 2,3. The on-shell version of Eq. (|TTD , 
appropriate to the choice /jS nf ^ = M h , is obtained by substituting l h —> L h = ln(M|/A 2 ) 
and c n Cn for N = 2,3. Analogously to the case of //w) = its inverse, which gives 
\( n f) in terms of A^ and Mh, then follows through the replacements A <->■ A'; (3^ P'n 
for N — 0,..., 3 ; and Cn —> — Cn for N = 2,3. 

In conclusion, we have extended the standard description of the strong coupling con- 






stant in the MS renormalization scheme to include four-loop evolution and three-loop 
matching at the quark-flavour thresholds. As a by-product of our analysis, we have set¬ 
tled a conflict in the literature regarding the two-loop matching conditions fllO], |TT|. These 


results will be indispensible in order to relate the QCD predictions for different observ¬ 
ables at next-to-next-to-next-to-leading order. Meaningful estimates of such corrections 
already exist 
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